Dissipative dynamics of vortex arrays in trapped Bose-condensed gases: neutron stars 
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We develop a theory of dissipative dynamics of large vortex arrays in trapped Bose-condensed gases. 
We show that in a static trap the interaction of the vortex array with thermal excitations leads to a 
non- exponential decay of the vortex structure, and the characteristic lifetime depends on the initial 
density of vortices. Drawing an analogy with physics of pulsar glitches, we propose an experiment 
which employs the decay induced heating of the thermal cloud as a tool for a non-destructive study 
of the vortex dynamics. 
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Superfluidity is one of the striking manifestations of 
quantum mechanics in interacting many body systems. 
The discovery of Bose-Einstein condensation in dilute 
atomic gases provides a novel example of super- 

fluid systems. Until recently, the evidence for superflu- 
idity in these systems has been fairly indirect and was 
based on either excitation spectrum or measurements of 
critical velocity (see for an overview of extensive ex- 
perimental and theoretical studies). A promising way of 
studying superfluidity in trapped gases is the observation 
of quantized vortices 

Aside from being a textbook example of persistent cur- 
rents, quantized vortices are believed to be widespread 
in Nature. They play an important role in breakup of 
superfluidity and superconductivity More compli- 
cated structures, such as arrays of vortex lines, appear 
in superfluids rotating with a sufhciently large angular 
velocity. This may occur in superfluid cores of rotating 
neutron stars, where the dynamics of vortices is thought 
to explain irregular features in pulsars signals (so called 
glitches, see for review). 

Vortex arrays in trapped Bose-condensed gases have 
recently been observed at ENS . Vortices were formed 
in the course of Bose condensation by stirring the gas 
sample with laser light. Depending on the angular fre- 
quency of the stirring beam, it was possible to create con- 
figurations with different number of vortices which form 
crystal-like structures. The vortex array is visualized by 
a destructive measurement of the gas cloud after ballistic 
expansion. The experiment shows that after the stirring 
ends, the vortex configuration disintegrates in a clearly 
non-exponential fashion. 

The idea to model superfluid interiors of neutron stars 
in experiments with vortices in liquid helium has a long 
history (see for an overview). Following this anal- 
ogy, vortex experiments in dilute gases may serve as an 
"Earth-scale" lab for studying the internal dynamics of 
rotating neutron stars. Namely, the Bose condensate 
filled with vortices can be used as a model for the su- 
perfluid core of the star. The thermal cloud interacting 
with the condensate can play a role of the non-superfluid 
crust of the star. We develop a theory explaining how 
the interaction between the two components of the sys- 
tem makes the vortices move towards the border of the 



condensate, where they decay by creating a sudden burst 
of thermal excitations (a "glitch"). We propose an exper- 
iment which employs the heating of the thermal cloud as 
a tool for a non-destructive study of the vortex dynamics. 

For simplicity, we consider a large array of vortex lines, 
emerged in a long cylindrical magnetic trap characterized 
by the radial trap frequency lu. Further on we assume 
that the angular velocity of the trap rotation 3> ilc, 
where flc is the critical velocity corresponding to the ap- 
pearance of the first vortex. In this case the equilibrium 
number of vortices N ^ 1. We assume that the stirring 
beam is on for a sufficiently long time so that the system 
reaches thermal equilibrium in the rotating frame. The 
equilibrium configuration of the system in the rotating 
frame can be obtained by variational analysis of the free 
energy functional 
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where V{r) = is the trapping potential in the 

radial direction, g — -iTTh^a/m, m is the the atom mass, 
fi is the chemical potential, and no(r), and 



M 



dr mnor x 



are, respectively, the superfiuid density, velocity and 
the angular momentum of the condensate in the labora- 
tory reference frame. The expression for the free energy 
given by Eq. (jl|) implies that the condensate is sufficiently 
dense: the mean interpaticle interaction fj. is much larger 
than the level spacing in the trap hu (Thomas-Fermi 
limit). 

In the presence of a vortex array with the vortex den- 
sity n, the superfiuid velocity satisfies Maxwell-like equa- 
tions 101 
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where z is a unit vector along the axis of the cylinder , and 
where k — Ti/m is the vortex circulation. The superfiuid 
velocity inside a cylindrically symmetric vortex array can 
be found from Eq.(ph by using the Stokes theorem. It 
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has no radial component (i.e. it satisfies the necessary 
boundary condition at the border of the condensate) and 
is given by 



v,(r) = 27ri 



r' dr'n{t, r'), 



(3) 



where f is the unit vector in the direction of r. 

The equilibrium configuration of no(r) and mini- 
mizes the free energy J-. The variational solution has to 
be found within the class of functions Vg corresponding 
to a given total number of vortices. Variation of the free 
energy (|^) with respect to the vortex velocity shows that 
in thermal equilibrium the velocity of the superfluid flow 
imitates the flow of a normal liquid with the angular ve- 
locity il: Vs = [il X r]. The corresponding density of 
vortex lines is also spatially uniform and is given by H 
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Then, the minimum of the free energy is reached at 
the condensate density given by the following modified 
Thomas-Fermi expression: 



no(r) = no(0) ( 1 - 



where the transverse size of the condensate is given by 

2fi 



and fi = gno(O). We note that the angular velocity can 
not be very large: The maximum velocity of the super- 
fiow, which coincides with the velocity of the superfiuid 
on the border of the condensate is ^ ilRc and shou ld not 
exceed the velocity of Bogolyubov sound Cg ~ ^/Jl/m. 
We thus have 
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since otherwise the superfiow becomes unstable with re- 
spect to spontaneous creation of excitations (see Landau 
arguments in [^). This limits the maximum number of 
vortices inside the vortex array: N ^ Nmax — 
which is still a large number in the Thomas-Fermi limit. 
In what follows we will assume that the number of vor- 
tices is still suflticiently small 1 ^ iV ^ Nmax, so that the 
angular velocity of the condensate rotation is smal com- 
pared to the trap frequency. Therefore, we can neglect 
the size increase of of the condensate due to the rotation. 

When the rotation of the trap ceases (i.e. after the 
stirring light is switched off and the rotation of the ther- 
mal cloud stops, see ||lj]), the vortex array continues to 
rotate with the angular velocity O. The vortex config- 
uration is described by two continuous functions: the 
average vortex density n{t, r) and velocity Vy{t, r). Since 
the number of vortex lines is a locally conserved quantity, 
the introduced functions satisfy the continuity equation: 
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= 0. 
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With this "macroscopic" description of the vortex "mat- 
ter" , we take a qualitative view on the vortex dynamics, 
ignoring the issues of metastable vortex configurations 
and the nature of elementary excitations in the vortex 
array (see [|l2| for a more detailed discussion). 

The vortex velocity is related to the superfluid velocity 
through the Magnus law Q : 



Ps[v^ - Vs, k] = F, 
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where ps is the superfluid (mass) density and F is the 
friction force acting on the vortex line (per unit length). 
At finite temperatures, the force F originates from the 
scattering of thermal excitations from moving vortices 
and consists of two terms: 
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Here D and D' are the longitudinal and transverse fric- 
tion coefficients, respectively. Both of them are temper- 
ature dependent and D' = p„, where p„ is the density 
of the normal component. The discussion of the friction 
coefficient D can be found in ||lj,|l5|. The solution of 
Eq-(0), with Vg from Eq.(^, can be represented in the 
form: 



v„ = w^^^f + w'"^) [k X f], where 
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are the radial and tangential components of the veloc- 
ity field, and p = -I- p„ is the total density. Typi- 
cally 13 <C K/9 and p p^ (see below). Hence, we have 
vi^"^ « Dvs/np <C Vs and v\f ^ « Vs- This means that in 
the presence of a finite friction {D ^ 0) the vortices cir- 
culate around the trap center with velocity ■^s and slowly 
spiral out towards the walls of the vessel. At the bor- 
der of the trap the condensate density vanishes and the 
vortex decays, creating a burst of excitations and thus 
transferring its energy to the thermal cloud. 

Assuming that the vortex density does not depend on 
the polar angle, i.e. n(f,r) = n{t,r), and using Eqs.(||) 
and (^), we rewrite Eq.(|g) in the form 



dn{t, r) TTKa d 



dt 



r dr 



n(t,r) / r'dr'n{t,r') = 0, (10) 



where a = 2kPsD / {k^ p^ -\- D"^) ^ 2D/ np. Although this 
equation is not local, the density at the point r depends 
only on the density profile inside the cylinder of radius r. 

The solution of Eq.(lO), with the initial density pro- 
file given by Eq.(^, remains spatially uniform at any t 
{n{t, r) — n{t)) and satisfies the equation 



dn{t) 
dt 



+ TTKan^it) = 0. 



Given the initial condition n(0) — Ueq, the total number 
of vortices N{t) = ■nR'^n{t) decays non-exponentially: 
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Here N{Q) = irR^neq, and the characteristic relaxation 
time Tfl of the vortex structure depends on the initial 
density 
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This result makes sense even if the number of vortices is 
not large. For example, in the case of a single vortex, i.e. 
for = 1, we have n^q ^ 1/^^, where R is the size of the 
vessel. Then, by substituting this density into Eq.(p^) 
we find that ^ nah/mR'^ . This differs only by a 
logarithmic factor from the life-time of a single vortex 
in a static trap [|l5|. The similarity of the results is not 
surprising, since the scenarios of dissipative dynamics in 
both cases are qualitatively the same. 

To be more specific, we present explicit expressions 
for Tn in a weakly interacting Bose-condensed gas. First 
we consider the limit of very small temperatures T <C 
H, where fi = inh^ano/m is the chemical potential of 
the gas, no is the condensate density, and a is the 2- 
body scattering length. Then, since ps ~ mno, and D cx 

^^5/2j.5/^7/2;j3 ^^^^ D ^ Kp and 
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In the opposite limit (/i ^ T) the calculation of the 
friction coefficient gives D k. Q.Ikw^'/'^T p^/"^ /h^ [ p^ . 
Then, assuming p„ ^ p, we obtain 



1.6.5/2^(^^^3)1/2?:^ 

Tri?^ p 
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In a spatially homogeneous gas the condition p„ ^ p 
used to derive Eq.(p^ implies T <C Tc, where Tc is the 
temperature of Bosc-Einstein condensation. In a har- 
monic trap characterized by the frequency w, thermal 

particles and the condensate occupy different volumes: 
Vt = (2T/tow2)3/2 j^^j ^ (2^/mw2)3/2 ^ y^^ j.g_ 

spcctively. Hence, the density of the thermal compo- 
nent can be much smaller than no even at temperatures 
close to Tc- This condition requires the fraction of Bose- 
condensed particles, which is ~ (T — Tc)/Tc, to be larger 
than VcIVt = [p/Tcf^ « 1 (see §). 

Given the values of experimental parameters | [Tl| ] : uj = 
169 X 2'Ks~^ and the number of particles 10^, we find that 
p w 100n_ft' and R = (2/i/mw^)^/^ « bpm. At temper- 
atures T p the size of the thermal cloud matches the 
size of the condensate. This prevents further evapora- 
tive cooling and complicates the temperature measure- 
ments. For this reason, we estimate the relaxation time 
for T ^ p, where both Eqs. (|l|) and (|l|) give approxi- 
mately the same result: 



For the array of 6 vortices we estimate t/j « 0.7s, which 
is in good agreement with the measured value w Is. 

Vortices in trapped Bose-condensed gases are very 
small and can not be visualized non-destructively by op- 
tical means. This is similar to the situation in astro- 
physics, where indirect evidence for the vortex dynamics 
was obtained by observing irregular features in signals of 
pulsars associated with rotating neutron starts |10| . In- 
teriors of neutron stars may contain superfluid nuclear 
matter sustaining huge vortex arrays. Every time a vor- 
tex reaches the crust of the star, it decays and the sudden 
release of its energy is detected on Earth as a glitch in 
an otherwise periodic signal of the star. 

In trapped gases the vortex decay transfers vortex en- 
ergy to the thermal cloud and hence increases the gas 
temperature. As we have noted, the superfluid flow with 
a large number of vortices imitates the rotation of a nor- 
mal liquid. Therefore the energy of a vortex array can 
be considered as an energy of rotation of such a liquid: 
E = /il^/2, where / = TrmnoRzR'^ /2 is the momentum 
of inertia of a condensate, Rz = y^2p/mujz and ujz are 
the condensate size in axial direction and the axial trap 
frequency, correspondingly. Given the heat capacity of 
the gas cloud C oc /h^uj'^ujz, we find that the temper- 
ature increase after the decay the vortex array: 
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The quantity AT can be large, at least for the tempera- 
tures well below the Bose-Einstein condensation temper- 
ature Tc- The heating of a Bose-condensed gas, associ- 
ated with the decay of the vortices, can be detected by 
a non-destructive technique. Then, for instance, by plot- 
ting the temperature of the thermal cloud as a function 
of time, one can confirm the decay scenario described by 
Eq. (ITTl) . This may provide a new tool for the study of 
the vortex dynamics and, at the same time, serve as a 
laboratory model of a rotating neutron star. 

In conclusion, we have developed a theory of dissipative 
dynamics of large vortex arrays. The calculated charac- 
teristic lifetime of the vortex structure in a static trap 
depends on the initial density of vortices and agrees well 
with the recent experiment |ll|| . Putting forward the 
analogy between the dissipative dynamics of vortices in 
trapped gases and the physics of neutron stars, we pro- 
pose an experiment, where vortices can be indirectly and 
non-destructively studied by measuring the heating of 
the thermal component 
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